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Abstract 

We consider the performance of the ATLAS and CMS searches for events with missing transverse 
energy, which were originally motivated by super symmetry, in constraining extensions of the Stan- 
dard Model based on extra dimensions, in which the mass differences between recurrences at the 
same level are generically smaller than the mass hierarchies in typical super symmetric models. We 
consider first a toy model with pair-production of a single vector-like quark U\ decaying into a 
spin-zero stable particle A\ and jet, exploring the sensitivity of the CMS ay and ATLAS m e ff 
analysis to M\j x and the U\ — A\ mass difference. For this purpose we use versions of the Delphes 
generic detector simulation with CMS and ATLAS cards, which have been shown to reproduce the 
published results of CMS and ATLAS searches for super symmetry. We then explore the sensitivity 
of these searches to a specific model with two universal extra dimensions, whose signal is dominated 
by the pair production of quark recurrences, including searches with leptons. We find that the LHC 
searches have greater sensitivity to this more realistic model, due partly to the contributions of sig- 
natures with leptons, and partly to events with large missing transverse energy generated by the 
decays of higher- level Kaluza-Klein recurrences. We find that the CMS ay analysis with ~ 5/fb of 
data at 7 TeV excludes a recurrence scale of 600 GeV at a confidence level above 99%, increasing to 
99.9% when combined with the CMS single-lepton search, whereas a recurrence scale of 700 GeV 
is disfavoured at the 72% confidence level. 
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1 Introduction 

One of the most powerful tools in the LHC quest for physics beyond the Standard Model is 
the search for events with missing transverse energy (MET), such as might be carried away by 
invisible dark matter particles. The prototypes for models predicting MET events have been 
provided by supersymmetry with R-parity, and the published CMS and ATLAS searches 
for MET events have been designed largely with supersymmetric models in mind pp. They 
are particularly sensitive to models with a large hierarchy between the masses of strongly- 
interacting sparticles, which are the most copiously produced at the LHC, and the dark 
matter particle, which has often been assumed to be the lightest neutralino. However, the 
sensitivities of some of these searches to supersymmetric models with smaller mass hierarchies 
have also been explored [2H3]. 

There are, however, many other classes of extensions of the Standard Model that also 
predict MET events, some with different generic features from the supersymmetric models 
with a large hierarchy of masses for which the initial LHC MET searches were optimised. 
Examples are models with universal extra dimensions (UEDs) compactified in a space with 
characteristic scale R ~ 1/TeV [HE]- Since this compactification scale is much smaller 
than the typical unification scale in supersymmetric models, the masses of the Kaluza-Klein 
recurrences of known particles are renormalised by relatively small amounts, and the mass 
differences between, e.g., the recurrences of strongly- and weakly- interacting particles at the 
same level are much less than would be the mass differences between strongly- and weakly- 
interacting sparticles in models with high unification scales. If such an extra-dimensional 
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model has an exactly-conserved Kaluza-Klein K-parity [6j[7] , the lightest recurrence is stable, 
and hence present as a dark-matter relic from the Big Bang and liable to be produced in MET 
events at the LHC. A related class of extensions includes Little Higgs models (HIE] , which can 
be thought of as minimal four-dimensional versions of deconstructed extra dimensions |10j . 
The inclusion of a T-parity pX] can stabilise the lightest neutral partner of the gauge bosons 
and thus also offer a dark matter candidate. 

The sensitivity of the initial ATLAS m e ff search to supersymmetric models with differ- 
ent mass hierarchies was studied in [2], and the sensitivity to the soft QCD jet-matching 
procedure and the sensitivities of other searches by both ATLAS and CMS (ax, vriTii razor 
variables and monojet) have been studied in [3]. In both cases, simplified models containing 
only near-degenerate gluinos, squarks and an LSP have been used. However, there are im- 
portant differences between these and extra-dimensional scenarios. In particular, the spins 
of the Kaluza-Klein recurrences are different from those of the hypothetical supersymmetric 
partners of Standard Model particles, affecting both production cross sections and decay pat- 
terns. For example, the decay chains necessarily include a sizeable fraction of leptons from 
the decays of the recurrences of the SU(2) gauge bosons, so searches containing leptons and 
MET must also be considered [T2|[T3] (this is also true for supersymmetric spectra contain- 
ing light sleptons). More importantly, the production and decays of higher-level recurrences, 
which are even under the Dark Matter parity, may also make important contributions to 
prospective MET signatures. 

In this paper we study various aspects of these potential features in searches for Kaluza- 
Klein excitations via MET searches at the LHC. Our first step is to define a minimal sim- 
plified toy model comprising a single excitation of the u quark, U\, that decays exclusively 
Ui — > u + Ai, where Ai is a neutral, weakly-interacting scalar particle. This toy model 
should encompass important features of the signal in realistic UED and Little Higgs models, 
which is dominated by the pair production of the lightest K- or T-odd recurrences of the 
light quarks. Here we consider a scalar A\, which is typical of 6-dimensional models, whereas 
a vector A\ would be expected in 5-dimensional and Little Higgs models. The toy model is 
also closer to the simplified models used to study compressed supersymmetric spectra, thus 
allowing us to make a direct comparison of the acceptances of the searches on the two kinds 
of signals 0- We study the prospective sensitivities of the searches looking for MET and jets 
in the (M^, M4J plane of this toy model, using variants of the Delphes generic detector 
simulation with CMS and ATLAS detector cards that have been validated for MET searches 
in the contexts of supersymmetric models. In particular, we use simulated signal events to 
analyse the efficiencies of the CMS ax, monojet, single-lepton (Lp), opposite-sign (OS) and 
same-sign (SS) dilepton analyses, and of the ATLAS m e ff search using the 2011 data at 7 
TeV, which have now been finalised. For definiteness, and because it has greater sensitivity 
than the ATLAS m e ff search and monojet searches, we focus on the CMS «t analysis, with 
a view to understanding how to optimise the strategies for future searches for Kaluza-Klein 
excitations. 

The main goal of this preliminary analysis is to validate our implementation of the exper- 

1 A similar toy model has been used in [14] to study the reach in Little Higgs models with T-parity: in 
our study, however, we focus on a region of parameter space with smaller mass splitting. 
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imental cuts, as a step to extending the analysis to a complete and realistic UED-type model. 
Several possibilities for such models have been proposed in the literature: in the present work 
we focus on the simplest possible model with a pair of universal extra dimensions, which 
is the minimal model with an exact K-parity proposed in [15]. The main advantage of this 
model is that it has an exact symmetry, which is a property of the geometry of the compact 
space and does not require any additional assumption on the interactions localised at the 
singular points of the space. We chose this model also for the following two reasons: the 
loop-induced mass splitting in each tier of recurrences is the smallest among the proposed 
models, thus offering an ideal arena to probe the sensitivity to compressed spectra. Also, 
the preferred mass range indicated by the Dark Matter abundance is 700 < ttikk < 1000 
GeV [IE] , which implies large production rates at the LHC and the possibility of probing the 
preferred masses with the present data. In contrast, UED models in five dimensions prefer 
larger values of m KK ~ 1.5 TeV [TJ. 

We have implemented the model in MadGraph5 [IT], and performed an analysis of its 
parameter space, including the full pattern of decays of the lowest even and odd tiers of 
excitations. The production cross section is dominated by the pair production of the lightest 
K-odd recursions of the light quarks, which roughly correspond to the toy model, while other 
contributions, specific to extra-dimensional models, come from the K-odd scalar gluon and 
the higher K-even recurrences. We use the complete simulation to check first how closely 
the events from the K-odd quark pair production are described by the toy model. We then 
compare the full set of channels, and study in detail the effect of all the searches, including 
MET searches involving leptons, in order to identify the production channels and searches 
that provide the main contributions to the combined bound. This exercise makes it possible 
to compare the bound with those obtained from other searches without MET, primarily the 
search for resonant dileptons coming from the decay of even recurrences of the neutral gauge 
bosons [IS]. It also provides some indications how to modify the supersymmetry searches to 
increase the sensitivity to UED-type models. In this paper we only include the final searches 
based on the 2011 dataset at 7 TeV, and we postpone the analysis of the preliminary 2012 
results and a study of the performance of alternative searches in a follow-up publication. We 
expect the inclusion of the 8 TeV data to certainly improve on the bounds, however without 
changing significantly the conclusions of this paper. 

The layout of the paper is as follows. In Section 2 we introduce the toy model with a single 
u excitation U\ and a stable neutral scalar A±, describe our implementations of the CMS 
ax and ATLAS m e ff MET analyses using Delphes, and go on to discuss the sensitivities of 
these searches in the (M^, M^J plane that characterises the toy model. Then, in Section 3 
we review the simplest model with two UEDs, and discuss the sensitivity of the CMS MET 
searches, including single- and dilepton searches, to this more complete scenario. We find 
that the LHC searches are more sensitive to this more realistic model, thanks partly to 
the decays of higher-level Kaluza-Klein recurrences that generate events with large missing 
transverse energy, and partly to the contribution of the single-lepton signature. Based on 
these results, in Section 4 we make suggestions how the present CMS and ATLAS MET 
search strategies could be modified to be more sensitive to such a UED scenario, and in 
Section 5 we present a summary and conclusions. 
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2 Getting Started 



2.1 The Toy Model 

We describe now the toy model we use to exemplify the features we seek to explore in this 
paper. To keep things as minimal as possible, we consider only the recurrences of the light 
quarks: in generic UED models, they are nearly degenerate because their universal mass is 
fixed by the size of the extra compact dimensions, with only small splittings generated by 
loop corrections. In Little Higgs models, the masses are arbitrary, but degenerate spectra 
are preferred in order to minimise corrections to flavour observables. 

The most minimal realisation of the toy model consists of only two new particles: a 
heavy up-type quark, U±, and a neutral stable scalar, A±. In the following, we assume that 
Mjj 1 > Ma 1 i and that the only allowed decay channel for the new quark is JJ\ — >■ uA\ with 
coupling cjj 1 a 1 u = 2e/3c^. The couplings has been chosen to match the coupling of the 
hypercharge gauge boson, as one usually obtains in realistic models. However, the precise 
value of the coupling is not relevant for the study, as long as the decay is prompt and the 
width small enough compared to the experimental accuracy on the MET and jet energy 
measurements. The neutral scalar cannot decay into any lighter state, and is therefore a 
Dark Matter candidate, invisible to the detector. The only free parameters of the toy model 
are the heavy quark mass Mu 1 and the mass splitting AM = Mjj 1 — Ma 1 ■ In the following 
we will be interested in near-degenerate spectra, and we consider mass splitting below 100 
GeV0. 

2.2 Simulation Details 

The processes for pair production of new particles in the toy model can be grouped into 
three classes: 

( QiQi \ 

PP — ► { Qi-A-i > + n jets where Qi — U\,Ui, (1) 
{ A X A X J 

with subsequent decay of the heavy quark U\ into the invisible scalar A\ and a quark u. 
The QQ class is by far the dominant one, because it includes the QCD production of U^Ui, 
which depends only on the mass of the heavy quark. Another relevant production channel 
in this class is U\U\, mediated by t-channel A\ exchange: this process may give a sizeable 
contribution at large M Vl masses, because the suppression due to the electroweak couplings 
would be compensated by an enhancement of the up-quark PDFs. However, this process 
also depends on the model-dependent UiA^u coupling: for this reason, in the following we 
consider only the QCD pair production and neglect all the electroweak processes. The events 
we are interested in are triggered by the presence of high-p^ jets, which can be originated 
either by the decays of the heavy quark U\, when the mass splitting is large or the heavy 

2 This toy model can be considered as a UED equivalent of the simplified supcrsymmctric model that only 
contains degenerate squarks and an LSP, with gluinos and non-coloured sparticles assumed to be decoupled. 
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Mass (GeV) 


250 


300 


350 


400 


450 


500 


550 


600 


Cross-section (pb) 


22.1 


7.83 


3.14 


1.38 


0.65 


0.324 


0.168 


0.091 



Table 1: Cross sections for U\U\ production at 7 TeV calculated in the toy model at NLO- 

nnll jzgnzy. 

quark is produced with a significant boost, or by additional QCD jets due to initial- or 
final-state radiation (ISR or FSR). For the latter reason, we need to incorporate events 
with additional hard QCD jets and carefully take into account the matching of the hard jet 
generation with the soft QCD jet incorporation. 

The toy model signal has been simulated with MadGraph5 [T7], interfaced with PYTHIA [21] 
for parton showering and hadronisation, and Delphes [22] for the detector simulation of both 
ATLAS and CMS. The simulation has been conducted scanning over Mu t and AM in the 
following ranges: 

M Vx = {250,600} GeV in steps of 50 GeV ; 
AM = 5,20,35,50,75,100 GeV. 

The cross sections calculated at NLO-NNLL for UiUi production at 7 TeV for the different 
values of My 1 are shown in Table [U We generated 300,000 events for each of these choices of 
Mu l and AM to ensure our equivalent integrated luminosity was larger than that collected 
at 7 TeV by the ATLAS and CMS experiments. The region with smaller mass splitting has 
been explored in more detail because, in generic UED models, particles in the same tier are 
almost degenerate and the small mass splittings arise only at loop level. 

Initial and final-state radiation (ISR and FSR) have been included both at the level of 
matrix-element generation and the parton shower, and the matching has been implemented 
using the "shower-fc^" scheme [23] encoded in Madgraph5: with this choice, the number 
of rejected events from Madgraph is low, enhancing the number of available events. The 
matching parameters xqcut (equal to qcut in the "shower-Zc^" scheme) has been tuned for 
every M^ — AM configuration in order to optimise the matching procedure for each process. 
Kinematical cuts on jets and missing energy have been set to minimal values at the generator 
level; tighter cuts have been applied after the detector simulation, at the level of the analysis, 
following the cuts used by the ATLAS and CMS collaborations. More details are provided 
in the next sections. 

2.3 Implementations of Delphes and Kinematic Distributions in 
the Toy Model 

We use for our analysis versions of the Delphes generic detector simulation package [22] 
with 'cards' that emulate the performances of the ATLAS and CMS detectors. Our imple- 
mentations of Delphes had previously been validated in analyses of supersymmetry searches 
using MET signatures [21]. The first step in that validation was to use SoftSUSY [25] to 
generate the spectra for various CMSSM points, then PYTHIA [2T| to generate event samples 
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for each of these points, followed by Delphes with the ATLAS or CMS detector 'card'. We 
then implemented the experimental cuts obtained from [I], and calculated corresponding 
confidence levels for these points. As shown in [24] for the case of the ATLAS MET analysis 
with 5 fb" 1 of data at 7 TeV [26], this procedure reproduced quite accurately the 95% CL 
bound quoted by ATLAS for the CMSSM, and we have performed a similar validation for 
the CMS MET search. The success of these Delphes validations in another extension of the 
Standard Model with a complex spectrum of massive states reinforces the expectation that 
Delphes simulations are suitable for models whose kinematic distributions lie well within 
the analysis search regions, such as for the extra-dimensional models considered here. We 
do not make any attempt to simulate SM background predictions, for which this is not the 
case. 

The ATLAS search [26] consists of 11 overlapping signal regions, which differ in the 
numbers of jets required and in the cuts on m e ff that were applied. In all cases, events 
with an electron (muon) with pj- > 20 GeV (10 GeV) are vetoed. We have implemented 
all these 11 regions and retained for each choice of the toy model parameters the region 
that sets the most stringent bound. The 11 regions have very different sensitivities to the 
toy model, and the two regions with the highest sensitivities (largest signal yields) are in 
general A' and C(900). Both of these require missing transverse energy E™ tss > 160 GeV, 
two jets with pr{juh) > 130,60 GeV respectively and E™ lss /m e ff > 0.25. In addition, A' 
requires m e ff > 1200 GeV whereas C(900) requires two more jets with prC/3, j&) > 60 GeV 
and m e ff > 900 GeV. Fig. CD displays various kinematical distributions (normalised to unity) 
for the ATLAS 5/fb analysis for fixed quark mass Mjj 1 = 600 GeV and varying the splitting 
Mjj 1 — Ma x from 100 GeV (blue) to 5 GeV (red). We see in the panels significant variations 
in the Nj ets , and hence and E^ tss , distributions as Mjj 1 —M^ varies. For larger splitting, 
the average number of jets reconstructed above the respective pt threshold increases due to 
the fact that the jets from the decay have larger whereas, for very small splitting, the 
high-pr jets originate mainly from QCD radiation. Also shown in the top right panel is the 
m e ff distribution after the ATLAS event selection for region A', again normalised to unity. 
We observe that the m e ff distribution is rather insensitive to the mass splitting, and similar 
conclusions apply to other values of the parameters of the toy model. 

The CMS collaboration has implemented various searches based on different kinematical 
variables: rriT2 [27], razor variables [28] and «r [29], which have similar sensitivities for the 
CMSSM. Here we focus on the 7 TeV «t based search using 5 fb _1 , because its documentation 
provides the details needed for our analysis. The kinematic variable oct, defined as 

M j3 jH* - MH* 

for a di-jet event, where j% is the less energetic jet, is very effective for disentangling QCD 
events, which typically have ax < 0.5, from events with large MET. Fig. [2] displays various 
kinematic distributions for the CMS olt analysis with fixed Mjj 1 = 600 GeV and varying 
Mu x — M J 4 1 . The distribution exhibiting the largest variations between the different values 
of Mu 1 — Ma x is that in missing Ht (MHt), whose difference is reflected in the MHt/Ht 
distribution (lower left) and in the olt variable itself (upper left). On the other hand, the 
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Ef 5 [GeV] H T [GeV] 



Figure 1: Various kinematical distributions (Nj ets , m e g, E™ 1SS and Ht) for the ATLAS 5/fb 
analysis for Mu 1 = 600 GeV and varying Mjj 1 — Ma 1 from 100 GeV (blue) to 5 GeV (red). 
All the distributions are normalised to unity. 
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distributions in Nj ets and Ht are more stable, though we do observe a slightly larger average 
number of jets for larger splitting, due to the contribution of the jets from the decay. 

2.4 LHC MET Constraints on the Toy Model 

We now apply the LHC constraints to the toy model, so as to gain insight into the typical 
limits one may obtain in a generic extra-dimensional model. Later we also discuss the 
differences with a specific realisation in a complete model including not only a heavy quark 
and a stable scalar, but the full set of Kaluza-Klein excitations, including not only the lowest 
tier but also higher-mass recurrences, which have different kinematical behaviours in general. 
We note that the signatures of this toy model are similar to a simplified model [30] used for 
supersymmetric searches [31], which contains squarks and the neutralino LSP and produces 
events with pair production of squarks followed by the decay q — » qxo- The two models 
differ in the spins of the new states, which alter the production cross section and influence 
the kinematic distributions of the decay products. 

As already remarked, in the toy model the production cross section is largely dominated 
by the QCD process pp —> U-JJi, which depends only on the mass of the heavy quark U\. 
Subleading contributions are given by processes like uu — > U-JJi, which are mediated by 
t-channel exchange of the stable boson A±. This is of electroweak strength, and has been 
neglected for this analysis. Although we estimate the experimental sensitivities and bounds 
with a single heavy quark, the results can be generalised approximately to any number of 
near-degenerate quarks, simply by multiplying the production cross section by the number of 
heavy quarks. The cross sections for quark-antiquark pair-production processes receive large 
corrections from QCD contributions at higher orders. These corrections must be taken into 
account when extracting reliable bounds on the masses of new quarks, and for this purpose 
we have computed cross sections at NLO+NNLL using the tool of Ref. (T9j[20] , with the 
results shown in Table [TJ 

In the left panel of Fig. [3] we display the efficiency found in our Delphes simulation of 
the CMS «t analysis for the toy model, as a function of the heavy quark mass Mjj 1 and the 
mass difference Mjj x — Ma± ■ We see that the experimental efficiency for the toy model is of 
the order 10~ 3 to 10~ 2 , which is similar to the efficiencies for supersymmetric models with 
compressed spectra found in [2]. Generally speaking, we see that the efficiency is reduced 
for smaller mass differences, essentially because these lead to smaller amounts of MET for 
any fixed value of Mu l , and is increased for larger values of Mjj 1 , which tend to yield more 
MET for any fixed value of M Vl - M Al § 

The number of signal events expected to be found in the CMS «t search is the product 
of the efficiency factor shown in Fig. [3] with the cross section shown in Table [TJ for a given 
integrated luminosity. However, evaluating the expected significance of any signal requires 
a comparison with the background in each bin and its expected level of fluctuation, which 
we take from the CMS analysis. In the right panel of Fig. [3] we display the confidence level 
we calculate using the CL S statistic for exclusion by the CMS a? analysis in the toy model, 

3 Typical differences between the masses of the quark recurrences and the LKP in the realistic model 
studied later are ~ 20 to 40 GeV, well within the range displayed in Fig. [3] see Figf4] 
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Figure 2: Various kinematical distributions (olt, Nj e t 8! MHt/Ht and Ht) for the CMS ar 
analysis for Mu 1 = 600 GeV and varying M\j x — Ma x from 100 GeV (blue) to 5 GeV (red). 
All the distributions are normalised to unity. 
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Figure 3: Efficiency (left) and exclusion confidence level (right) of the CMS ax analysis for 
the toy model, as a function of the heavy quark mass Mjj l and the mass difference Mjj 1 —M^ ■ 



as a function of the heavy quark mass Mu 1 and the mass difference Mu x — ■ Generally 
speaking, the upper limit weakens as the mass difference M Vl — M& 1 decreases, and also as 
Mu x increases. The weakening of the limit as Mjj l — decreases reflects a lowering of the 
selection efficiency, as seen in Fig. El whereas the increase in the efficiency as Mu 1 increases 
does not counteract the rapid decrease in the cross section. We see that Mu x = 300 GeV 
is excluded beyond the 99% CL for any value of Mjj 1 — Ma ± , and that Mjj 1 = 350 GeV is 
excluded at the 95% CL for mass splittings above and including — Ma 1 > 50 GeV. A 
similar exercise for the ATLAS experiment yielded slightly weaker results, mainly due to the 
lower selection efficiency for toy model events, so we concentrate on the CMS ay search in 
the following 0. 



3 A More Realistic Model 

In this Section we introduce and analyse a more realistic model with Universal Extra Dimen- 
sions (UED), based on a two-dimensional orbifold of the Real Projective Plane (RP 2 ) |15j . 
There are two main reasons for the choice of this particular space. In the first place, it is the 
only orbifold in 5 or 6 flat dimensions where an exact Kaluza-Klein K-parity can be defined. 
In other realistic models proposed in the literature [S1E3], generic interactions localised on 
the fixed or singular points of the orbifold will in general break the Kaluza-Klein parity of 
the model explicitly. Such localised interactions correspond to higher-order operators, so 
their absence would require special assumptions about the ultra-violet completion of the 
extra-dimensional model. On the other hand, on the RP 2 the Kaluza-Klein parity is exact, 

4 We also found that monojet searches (32] gave significantly weaker levels of exclusion, so we also do not 
discuss them further. 
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whatever the properties of the ultra-violet physics. In addition, the mass splitting between 
different states in the same tier, which are generated mainly by loop corrections, turn out to 
be rather small compared to other UED models proposed in the literature [34, 35J . Therefore, 
this scenario is a natural representative of the class of models we want to explore in this 
paper. The second reason why we focus on this model comes from the phenomenology of 
the Dark Matter candidate: in order for the relic abundance to match the determination by 
WMAP and other measurements, the mass scale of the new states is preferred to lie in the 
range 700 < uikk < 1000 GeV [16]. Higher values of the masses would imply an overabun- 
dance of matter in the Universe, and are therefore excluded by cosmology 0. This mass range 
is ideal for LHC searches, as it offers production cross sections large enough to be probed with 
the present amount of data. Other models in the literature prefer very different mass ranges: 
in the minimal 5-dimensional model, the preferred range is 1250 < m^K < 1500 GeV [36J, 
thus too high to be accessible at present, whereas for the other 6-dimensional model based 
on the chiral square, the preferred range is very low, 180 < ttikk < 200 GeV [37], thus likely 
to be excluded already. 



3.1 A Minimal Scenario with Two Universal Extra Dimensions 

The RP 2 orbifold is defined via a discrete symmetry group of a two-dimensional Euclidean 
spaced describing the two extra space co-ordinates x 5 ,x e . It consists of a reflection, r, and 
a glide, g, whose actions on the co-ordinates are: 

r({x 5 ,x 6 }) = {-x 5 ,-2 6 }, (3) 
g({x 5 , x 6 }) = {x b + tiR 5 , -x 6 + tiRq} . (4) 

We note that g' = g * r defines a second glide that flips the sign of £5, and that, squaring 
the two glides, one obtains translation symmetries along the two co-ordinate directions: 

9 *g({x 5 ,x 6 }) = {x 5 + 27rR 5 ,x 6 } , (5) 
g' * g'{{x 5 ,x 6 }) = {x 5 ,x 6 + 2nR 6 } . (6) 

The fundamental space of RP 2 can therefore be thought of as a subspace of an orthogonal 
torus with radii R5 and Rq. 

The spectrum of the model is therefore the same as a torus, up to a projection defined 
by the orbifold symmetries: each tier of resonances can be labelled by two integer numbers 
(n, m) that represent units of the quantised momenta along the two extra directions. The 
bulk mass is therefore given by 

2 2 

2 n m (~\ 

m (n,m) ~ + ^2 ' 

5 We note, however, that by tuning the mass of the Higgs recurrences, it is possible to find a "funnel" 
region, where the annihilation cross section is dominated by resonant s-channel heavy Higgses, as in the case 
of supersymmetry, and the mass scale can be pushed up to ~ 1.5 TeV, at the price of fine tuning in the 
parameter space. 

6 The RP 2 orbifold can also be defined as a sphere with the identification of antipodal points. A UED 
model based on this space has been proposed in [35], and has very different phenomenology. 
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where the two inverse radii define the mass scales in the model. The field content of each 
tier of recurrences depends on the field content and on the parity assignment on each six- 
dimensional field. 

In the following, we focus on a UED model in which a bulk field is associated to each 
Standard Model (SM) field [15]. Note that chiral fermions are independent fields in the SM, 
therefore two 6-dimensional fields are associated to each SM fermion. Parities are assigned 
to each field in such a way that the zero-mode spectrum coincides with that in the SM. The 
conserved symmetry that is responsible for the stability of the Dark Matter candidate can 
be defined in terms of a translation 

Pk({x5,x 6 }) = {x 5 + irR 5 ,x 6 + iiR 6 } , (8) 

and the tiers have parity (— l) n+m under this symmetry. This is the only exact symmetry 
that is respected by all interactions in the model [18]. The degeneracy within the tiers of 
recurrences with the same values of (n, m) is lifted by loop corrections and by the contribution 
of the Higgs VEV. 

In the following we are mainly interested in the lightest tiers: the odd tiers (1,0) and 
(0,1), and the even tiers (2,0) and (0,2). Each odd tier contains a gauge scalar boson 
corresponding to each SM gauge vector boson, and a vector-like fermion for each chiral 
Standard Model fermion. To simplify the notation, we label the odd states with a subscript 
"1" : the complete field content of the odd tiers is therefore 

gauge scalars =>■ A\ y Z\, W\ , G\ ; 

fermions q is , q 1D , lf L , lf R , v\ , 

where q — u,d, s, c, b, t labels the quarks, I — e,/i,r the charged leptons. The odd states in 
these tiers can only decay to the lightest state in the tier, which is the Dark Matter candidate 
A±, plus a SM particle. The even tiers contain a vector gauge boson corresponding to each 
Standard Model gauge boson, a vector-like fermion for each Standard Model chiral fermion 
and a massive scalar Higgs doublet. We denote even states with a subscript "2": 

gauge vectors =>■ A^, W^, G^; 

fermions q 2S , q2D, lf L i l tw v 2\ 
Higgs scalars =>■ H 2 , S^, S^ 1 . 

Here the are the recurrences of the Goldstone bosons of the SM Higgs, i.e., the charged 
component of the doublet and the imaginary part of the neutral field (pseudoscalar) . De- 
pending on the mass corrections, a state in an even tier may decay into a pair of odd (1) 
states, into an even state (2) plus a Standard Model particle, or into a pair of Standard 
Model particles. The latter would arise from a loop-induced coupling [39], but the rate is 
comparable to the other two processes, which are suppressed by the phase space available af- 
ter the loop-induced mass corrections. The MET signal comes typically from the production 
of odd states either directly or via the decays of an even state []. 

7 We note that there is another even tier that is potentially lighter that the (2,0) and (0,2) one, namely 
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The spectra of the four tiers we are interested in depend on two free parameters, namely 
the two radii -R5 and Rq. The loop corrections are also logarithmically sensitive to a cut-off, 
i.e., the scale where the 6-dimensional theory is no longer under perturbative control. In 
the following, we fix this scale to be 10 times the mass scale of the Kaluza-Klein states, 
which is the largest value allowed by naive dimensional analysis. Lowering the value of 
the cut-off would decrease the mass splitting, but also increase the sensitivity of direct 
Dark Matter detection experiments to the level that values of the cut-off below Am^K ar e 
already excluded [16]. The spectrum also depends on the higher-order operators, localised 
on the two singular points, which act as counter-terms for logarithmic divergences in the 
loop corrections. In the spirit of effective field theories, we assume that the contributions 
of the higher-order operators are negligible compared to the loop corrections. When writing 
the loop correction to the mass as m = m( n m )(l + 5), it should be noted that 5 has a mild 
dependence on the ratio R$/Rq. This is due to the fact that the main contribution to the 
loop correction is proportional to the Kaluza-Klein bulk mass of the state [IB] . The localised 
operators will modify the spectrum also by introducing mixing between the two odd and 
the two even tiers. However, while the even tiers can also mix via loop contributions, the 
odd ones only mix via contributions of the ultraviolet physics. In our phenomenological 
analysis, therefore, we need only consider mixing in the even sector, which can be rather 
large if the values of the two radii are close. This fact allows us to distinguish two separate 
regimes: one where the radii are almost degenerate, their values being similar within 10%, 
and one where there is a hierarchy between the two radii. The main difference between the 
two limits is in the preferred mass range [16]: in the degenerate limit, the mass scale tends 
to be smaller and is already excluded for all values of the cut-off by direct Dark Matter 
detection experiments. On the other hand, in the hierarchical case, only one pair of sets of 
states associated with the larger radius, say (1, 0) and (2, 0), are relevant, and the preferred 
mass range is 700 < m KK < 1000 GeV, with values of the cut-off below 4 excluded by direct 
Dark Matter detection experiments. 

In the following, we only consider the hierarchical limit: in this case, the phenomenology 
is dominated by the lightest tiers, which get their masses from the larger radius. The results 
therefore only depend on one parameter, m^K — min(l/ 'Rs,l/ 'Rq), which sets the overall 
mass scale. In Figure H] we show the mass corrections 5mm = m (i) ~ ^kk for the lightest 
odd tier as a function of m^K at one loop. The figure shows that the mass difference between 
each KK recurrence and the lightest one over the relevant mass range amounts to a few tens 
of GeV for coloured states (quarks and gluon recurrences), i.e., within the range studied 
previously in the toy model, and few GeV for leptons. In the right panel of the same figure 
we show the same mass corrections for the even tier, Sm^) — m (2) ~ ^ m KK- 

the (1, 1) tier. However, it can decay into a pair of Standard Model states only via localised higher-order 
operators, whilst the tree-level process (1, 1) — > (1,0) + (0, 1) is below threshold. We will not consider this 
tier any further, because it does not lead to MET signals. Some aspects of its phenomenology have been 
studied in [40]. 
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Figure 4: Mass correction for the lightest odd (1) and even (2) tiers of recurrences in the 
hierarchical limit. We show the correction with respect to the bulk mass, respectively tukk 
and IrriKK- The various colours label recurrences of: the electroweak gauge bosons in blue (A, 
W ± and Z , from bottom to top), leptons in green (singlets and doublets), light quarks in red 
(up and down singlets, doublet ant bottom- doublet) , top quarks in violet (singlet/ doublet), 
gluon in magenta and higgses in black (S2 , S2 and H 2 , where S are recurrences of the 
Goldstone bosons). 



3.2 Simulation of the RP 2 Scenario 

The MET signals we study in this work are generated by the decays of odd states into the 
dark matter candidate A\. In the following, we consider two classes of signals (^1,2 labels 
any particle in the odd or even tiers) 

1) pp X1X1, which is dominated by the production of quarks, and can therefore be 
compared with the toy model in the previous section; 

2) pp — y X2X2, where one of the even states decays into a pair of odd states and the 
other decays into energetic SM particles (jets, tops or leptons) thus providing the 
visible energy for the analysis. 

The processes in class 1) give the largest production rates, due to the smaller masses, 
and they have similar kinematics as the toy model. The processes in class 2) are particularly 
interesting, as they are absent in minimal supersymmetric and Little Higgs models, where all 
(or most of) the new states are odd under the dark matter parity. The complete RP 2 model, 
including the loop-induced vertices, has been implemented in MadGraph5 with the help of the 
FeynRules [41] package, and tested by comparing branching ratios and cross sections with 
the CalcHEP 3.0 [12] implementation. In the implementation we included the full one-loop 

8 Another potentially interesting class is pp — > X 2 — > X^Xi and pp — > X 2 — > xsmX 2 , where the resonant 
even state is a gauge vector (G2, W^, Z 2 or A 2 ) and its production takes place via a loop- induced vertex. 
We found that the matching procedure used in this work has difficulties with the latter process, where xsm 
is a jet, therefore including this process would entail calculations beyond the scope of this work. To be 
conservative, we did not consider single resonant production. 
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results for the masses of both even and odd recurrences, which depend only very mildly on 
the ratio of the two radii [18] , as noted previously The mass splittings between various states 
in the odd (1) and even (2) tiers are shown in Figure HI Since the mass differences between 
coloured and electroweak states are typically a few tens of GeV, production at the LHC is 
dominated by coloured states. The decay chains we are interested in are therefore initiated 
by heavy quarks and gluon recurrences. These chains can be very complicated, because they 
can involve the lighter W and Z recurrences and the leptons. For the even tiers, decays into 
a pair of resonant SM particles can also occur with a significant rate. Quarks, for instance, 
always decay into a SM quark plus an electroweak gauge boson in the same tier: this is true 
for both odd and even recurrences. The subsequent decays of the gauge bosons will involve 
further jets and, in some cases, leptons. 

The possible final states of the decay chains of the odd quarks, grouped according to the 
numbers of jets and leptons (electrons, muons and leptonically decaying taus), are pictured 
in Figure \5\ together with their relative rates, for a benchmark value m^K — 400 GeV. 
The relative rates change only slightly for larger masses. These plots show that the singlet 
recurrences qis decay predominantly into a single jet plus a Dark Matter candidate, and 
so correspond directly to the toy model. On the other hand, doublet recurrences qm yield 
significant rates for final states with leptons or multiple jets, so they reflect physics not 
included in the toy model. The top recurrences preferably decay into final states with 
leptons, however their production cross section is much smaller compared to the recurrences 
of light quarks, thus their impact in the phenomenology is limited. 

Similar plots for the quark even recurrences are shown in Figure [6j here we note large 
rates in final states without MET, due to the resonant decays of the heavy gauge bosons in 
the chains. We note also that here the MET is due to a pair of Dark Matter particles, so 
the kinematic distributions generated by the decays of these particles are very different from 
the toy model, as we see later. The branching ratios for the odd and even gluon recurrences 
are shown in Figure 

In producing signal events, we have considered only the most relevant production pro- 
cesses of heavy states, namely those due to QCD interactions. Their production channels 
can be grouped in three sets: 

GiGi GiQi QiQi (9) 

where i = 1, 2 and Qi stands for both quarks and anti-quarks. In addition to the usual QCD 
interactions, the model also contains processes mediated by the heavy gluon recurrence: for 
instance, significant contributions come from the process uu — > UiUi, mediated by t-channel 
exchange of the gluon G\. Such contributions are included in our simulation. We checked 
that the contribution of electroweak processes is negligible, amounting to a percent level of 
the QCD cross section at most. 

The framework for the simulation and jet matching is analogous to that used for the toy 
model. The main difference arises from the rich structure of chain decays of RP 2 states, which 
leads to a great variety of different final states. A simulation of production and decay of heavy 
RP 2 states within Madgraph would be computationally extremely demanding. Therefore, we 
have used an external tool to decay heavy particles in the events, taking into account the 
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Figure 5: Branching ratios for quarks of the (1,0) tier decaying to final states containing jets, 
leptons and missing energy. Only partners of the up and top quark are shown: the branching 
ratios of other quark recurrences are similar to those of the up recurrences. 
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Figure 6: Branching ratios for quarks of the (2,0) tier decaying to final states containing jets, 
leptons and missing energy. Only partners of the up and top quark are shown: the branching 
ratios of other quark recurrences are similar to those of the up recurrences. 
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Figure 7: Branching ratios for gluons of the (1,0) and (2,0) tiers. 
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Odd 


tier 






Even tier 






M KK (GeV) 


QxQx 


QiQi 


GiQi 


G\G\ 


Q2Q2 


Q2Q2 


G2Q2 


G2G2 


total 


400 


1,630 


7,440 


4,780 


418 


718 


159 


476 


43 


15,700 


600 


221 


531 


327 


18 


25.8 


2.6 


7.2 


0.4 


1,130 


700 


99 


179 


119 


5.7 


5.2 


0.36 


1.08 


0.05 


409 



Table 2: Production cross sections in fb at 7 TeV for the 8 classes discussed in the text, 
which are distinguished by the pairs of odd and even states produced. 

helicity amplitudes, namely the software BRIDGE [13], which decays particles in event files into 
every possible final state, taking into account spin correlations. This procedure is, however, 
an approximation because spin correlation effects between the Madgraph output and the 
BRIDGE input cannot be fully taken into account, but kinematical properties of final states 
are reproduced with a reasonable accuracy. In any case, the detector simulation smears 
physical observables and, moreover, we do not search for asymmetries or polarisation effects 
in the final states. To analyse a specific final state it is therefore enough to generate a large 
sample of events and select the subset of events only in the final phase of the analysis. 

In making our simulation we grouped similar production processes into 8 classes: GiG± 
and G 2 G 2 comprising pair production of gluons, G1Q1 and G2Q2 comprising the production 
of a gluon with a quark or anti-quark, and Q2Q2 comprising pair production of two 

quarks 0, and finally Q1Q1 and Q2Q2 comprising the production of a quark and antiquark. 
In all classes, we include all flavours of the heavy and SM quarks and, as in the toy model, 
processes with 1 and 2 additional jets are included. The matching procedure has been opti- 
mised for each class and each mass choice independently, following the procedure described 
in Section |2~2"1 We studied three benchmark points, with masses m^K — 400, 600 and 700 
GeV. The cross sections in each of the individual classes have been computed at leading order 
using our Madgraph implementation and are shown in Table [2j Unlike in the toy model, it 
was not possible to reliably estimate the QCD corrections to the cross section, due to the 
presence of many diagrams contributing to each production class, and to the presence of 
heavy gluons. Therefore, we conservatively used LO cross sections in this study. For each 
production class and mass point we generated 300,000 events, with a large majority passing 
the matching procedure with PYTHIA. 

3.3 MET Searches with Leptons 

For the analysis, we used the same validated Delphes implementation of the CMS and 
ATLAS detectors that we used to analyse the toy model. However, because of the presence 
of leptons in the final states, which were absent in the toy model, we now include, together 
with the MET analysis detailed in Section 12. 3[ typical supersymmetry-motivated searches 
based on leptons. Thus, following the philosophy of the previous section, we focused on 
the searches published by the CMS collaboration: together with the CMS ax analysis [2"§] . 
in this section we consider single- lepton [H], opposite-sign [15] and same-sign dilepton [15] 
searches. Similar searches have been performed by ATLAS, and their inclusion would not 

9 The production of two antiquarks is negligible. 
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0.00 
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0.4 
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0.00 


0.19 


0.00 


0.00 


G2Q2 


125.33 


1.57 


0.25 


36.74 


0.73 


0.15 


5.07 


0.11 


0.03 


2.39 


0.06 


0.01 


Q1Q1 


64.68 


8.51 


3.3 


2.17 


0.26 


0.11 


0.19 


0.04 


0.04 


0.19 


0.02 


0.01 


QiQ 1 


342.86 


24.06 


7.29 


12.53 


0.86 


0.24 


2.6 


0.29 


0.16 


0.5 


0.04 


0.03 


G1G1 


28.38 


1.7 


0.61 


0.65 


0.05 


0.02 


0.02 


0.00 


0.00 


0.01 


0.00 


0.00 


GxQx 


229.58 


20.78 


8.51 


6.29 


0.63 


0.27 


0.4 


0.08 


0.06 


0.24 


0.03 


0.02 


Total 


1085.67 


63.39 


21.24 


149.43 


6.31 


1.7 


22.23 


1.14 


0.45 


10.68 


0.46 


0.16 



Table 3: Numbers of events for each subprocess that pass the signal selections of the respective 
searches in 5/fb at 7 TeV. 



change the conclusion of our work. A study of the reach of the ATLAS multi-lepton searches 
for supersymmetric compressed spectra can be found in Ref. [17]. 

The single-lepton search relies on events with a single high-p-^ lepton together with jets 
and MET. Three different methods are used to discriminate the signal from the background, 
and the one we implement here is based on the "lepton projection" method that uses a 
variable L p |H] , which measures the component of the lepton transverse momentum that is 
parallel to that of the reconstructed W it originates from. For SM events, where the lepton 
comes from the decay of a W boson via a left-handed coupling, the lepton momentum tends 
to be aligned in the direction of motion of the W, resulting in a large value of L p . Therefore, 
a signal region is defined for L p < 0.15, while L p > 0.3 is used as a control region. Additional 
cuts on Ht and on the scalar sum of the transverse momentum of the lepton and transverse 
missing momentum, S l ^ pt , are also implemented (typically, Ht > 500 GeV). In the opposite- 
sign (OS) dilepton search [15], events with a pair of isolated leptons are selected, rejecting 
regions close to the Z peak in invariant mass. To suppress the ti background, 4 signal regions 
are defined with cuts on the missing transverse energy and Ht'- high-i?™^ (_E™ SS > 275 
GeV, H T > 300 GeV), high-# T (Eip iss > 200 GeV, H T > 600 GeV), tight signal region 
^ E miss > 275 GeV, H T > 600 GeV) and \ow-H T {EJ£ iss > 275 GeV, 125 < H T < 300 GeV). 
Finally, in the same-sign (SS) search [16], the presence of two isolated dileptons with the 
same charge is required, where the isolation is established by requiring that the scalar sum of 
the track momenta and calorimeter energies in a cone AR < 0.3 around the lepton track has 
to be less than 15% of the lepton momentum. Five signal regions are defined with various 
cuts on H T and EJp iss : Region 1 (H T > 80 GeV, EJp iss > 120 GeV), Region 2 (H T > 200 
GeV, EJp iss > 120 GeV), Region 3 (H T > 450 GeV, E™ iss > 50 GeV), Region 4 (H T > 450 
GeV, E™ iss > 120 GeV), and Region 5 {H T > 450 GeV, E™ iss > GeV). 

In Figure [S] we display the total number of events (Table ED that pass the selection cuts 
of the 4 searches in question for the 3 benchmark masses, showing in particular the relative 
contributions of the 8 classes of production subprocesses. We see that the numbers of events 
surviving the L p single-lepton cuts are typically ~ 10% of the numbers of events surviving 
the aT MET selection, with much smaller numbers for the opposite- and same-sign (OS and 
SS) dilepton searches. As a consequence, the latter contribute very little to the overall model 
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Figure 8: Bar-charts with the total numbers of events in 5/fb at 7 TeV that pass the sig- 
nal selections of the respective searches, showing the relative contributions of the individual 
classes of subprocesses. 
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Figure 9: The normalised (left) and MET distributions (right) in the toy model with 
M Vl = 400 GeV and AM = 20 GeV (red) and for the Q 1 Q 1 and Q 1 Q 1 subset of the full 
model with txikk = 400 GeV (black). 



sensitivity, which is dominated by the analysis with significant support only from the L p 
analysis. 

3.4 Contributions of Level (2,0) States to MET Searches 

We have already mentioned some reasons why the results in the realistic model should differ 
to some extent from those obtained in the toy model. For one thing, there are additional 
particles in the (1,0) tier, beyond the analogue of the U\ quark recurrence discussed in 
the toy model, leading to the more complicated decay chains with leptons discussed above. 
Secondly, in the realistic model, higher tiers are present: in particular, the even level (2,0) 
quarks play an important role in the full model. They contribute to subprocesses with 
different kinematics from the (1,0) tier quarks, giving distributions that differ from those of 
the toy model. 

In order to qualify and quantify the differences, in Figs. [9] and [10] we display two key 
kinematical variables used in the analysis, namely Hj> and MET, in the RP 2 model for 
m KK = 400 GeV, and in the toy model for M Vl = 400 GeV and AM = 20 GeV. As a first 
cross-check, Fig. [9] shows the sum of the QiQi and QiQi subprocesses (black) compared 
to the toy model histograms (red), with both distributions normalised to the same number 
of events. The similarity in the shapes are as expected, and the comparison validates the 
toy model approach. When compared to the full model in Fig. [10] however, we see a much 
stronger extension of the distribution (left) above 300 GeV with a notable bump around 
800 GeV, probably coming from resonant hadronic decays of the even (2) gauge bosons. 
We also see in the MET distribution (right) a peak at small values but also a stronger 
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Figure 10: The normalised Ht (left) and MET distributions (right) in the toy model with 
M Vl = 400 GeV and AM = 20 GeV (red) and in the full model with m KK = 400 GeV 
(black). Note the important contributions of the (2,0) states at larger Ht and MET, and the 
contributions of more complicated decay chains at smaller MET. 



extension above 200 GeV. These extensions of the Ht and MET distributions are due to the 
contributions of the second-tier (2,0) states, as is the low-MET peak. The extensions of the 
Ht and MET distributions imply that the sensitivities of typical supersymmetry-motivated 
MET searches at the LHC to the full model are greater than those to the toy model, as we 
now discuss. 

The impact of the (2,0) states on the searches is also shown in Fig. [Sj where we see that, 
in the case tukk — 400 GeV, the even states contribute ~ 40% of the events expected in the 
«t MET search and over 80% of the events expected in the L p single-lepton search. For the 
«t search, which contains the largest number of events, the importance of the even states 
decreases considerably for larger masses, dropping to about 10% and 5% for ttikk = 600 and 
700 GeV. The main reason behind is the decrease in the production cross sections of a pair 
of (2) states which is suppressed by a larger mass compared to the (1) states. In fact, the 
efficiencies of the ot search in the individual classes of production, listed in Table SJ show a 
mild decrease in all channels for increasing m^K- On the other hand, in the same Table H] 
we can see that, while the cross sections for pair production of the odd states at 700 GeV 
are between 1 5% of the values at 400 GeV, for the even states the decrease amounts to 
0.1 -T- 1%. The situation is different in the leptonic searches, like L p : in Fig. [8] we see that 
the even tiers are always very important, providing most of the signal events for all values 
of the masses. This effect can be understood by looking at the efficiencies of the L p search 
in Table HJ while the efficiencies decrease for larger masses in the classes of odd states, for 
even states the efficiencies increase at larger masses. This change in efficiencies suffices to 
compensate for the decrease in cross section. The increased efficiencies may be explained 
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Total 


cross section [fb] 


Q!t efficiency 




efficiency 


m KK [GeV] 


400 


600 


700 


400 


600 


700 


400 


600 


700 


Q1Q1 


7440 


531 


179 


0.92% 


0.91% 


0.81% 


0.034% 


0.032% 


0.027% 


G X Q X 


4780 


327 


119 


0.96% 


1.27% 


1.43% 


0.026% 


0.039% 


0.046% 


Q1Q1 


1630 


221 


99 


0.80% 


0.77% 


0.66% 


0.027% 


0.023% 


0.022% 


G\G\ 


418 


18 


5.7 


1.36% 


1.88% 


2.14% 


0.031% 


0.05% 


0.05% 


Q2Q2 


718 


25.8 


5.2 


6.76% 


4.76% 


4.7% 


2.05% 


2.63% 


3.35% 


G2Q2 


476 


7.2 


1.08 


5.27% 


4.36% 


4.6% 


1.55% 


2.0% 


2.7% 


Q2Q2 


159 


2.6 


0.36 


5.23% 


4.32% 


3.5% 


1.84% 


2.84% 


3.08% 


G2G2 


43 


0.4 


0.05 


4.8% 


4.0% 


4.1% 


1.27% 


1.5% 


2.0% 


total 


15700 


1130 


409 


1.4% 


1.1% 


1.0% 


0.19% 


0.11% 


0.08% 



Table 4: Cross sections and efficiencies of the ax and L p searches in the individual classes, 
and for the three recurrence scales considered, tukk — 400,600 and 700 GeV. 





m KK = 


400 GeV 


m KK = 600 GeV 


m KK = 700 GeV 




^total 


CL 


e total 


CL 


Ctotal CL 




1.4% 


100% 


1.1% 


99% 


1.0% 64% 


L p 


0.19% 


100% 


0.11% 


83% 


0.08% 38% 


OS 


0.03% 


87% 


0.02% 


3% 


0.02% 1% 


ss 


0.01% 


100% 


< 0.01% 


20% 


< 0.01% 5% 


Combination 


100% 


99.9% 


72% 



Table 5: The efficiencies etotai and the exclusion confidence levels CL for the signatures 
considered in the CMS 7 TeV data: the MET search, the L p single-lepton search, and the 
opposite- and same-sign (OS and SS) dilepton searches, as well as the combined CL. Results 
are given for the three recurrence scales considered, m^K — 400,600 and 700 GeV. 

by the presence of larger MET and leptons with higher p?, coming from the decays of the 
massive W% -^V^v. 

3.5 Constraints on the UED Scenario 

The results of the previous sections can be used to test the reach of the 7 TeV searches on 
the RP 2 UED model. Here we only use the results of CMS, in particular the ax [29], the 
L p single lepton [H], and same sign jl6] and opposite sign [15] (SS and OS) searches. The 
number of events that pass the selection cuts, listed in Table [3l can be combined with the 
experimental results to obtain a bound. We use the CL S statistic to calculate the exclusion 
confidence levels (CL) shown in Table [5]for ttikk = 400, 600 and 700 GeV, respectively. Also 
shown is the exclusion confidence level obtained by combining the different searches for each 
value of itikk- 

We see that the case m KK = 400 GeV is excluded with a very high degree of confidence 
independently by the «t, L p and SS analyses, and disfavoured to a lesser extent by the OS 
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analysis. This value of ttikk is therefore excluded very robustly. The case rriKK = 600 GeV 
is excluded at the 99% CL by the analysis alone, and at the 83% by the L p analysis, 
but not significantly by the dilepton analyses, which start losing reach at larger masses. 
Combining the analyses, we find that this value of tukk is also excluded robustly, at the 
99.9% CL. However, the case uikk = 700 GeV cannot be excluded: it is disfavoured at the 
64% CL by the ax analysis and at the 38% CL by the L p analysis, but the combined CL is 
only 72%, with negligible contributions from the OS and SS analyses. 

We can therefore conclude that the LHC MET searches (CMS) at 7 TeV require a bound 
on the ttikk mass in the RP 2 model between 600 and 700 GeV. It is interesting to compare 
this number with the expected bound from searches without MET: the most promising one 
is the bound on resonant di-leptons originating from the decays of A 2 and Z 2 , which are 
produced in the decay chain of even coloured states. Even though the events contain other 
particles, and the leptons are coming from two resonances, the number of events can be 
directly compared to the exclusion for a single Z 1 — > l + l~ . It was shown in [TB] that the 
corresponding bound is m KK > 575 GeV, and other channels such as W% — > l^v and jet 
resonances are less competitive. This comparison shows that searches with MET are already 
competitive with searches that do not require the presence of MET, and seem already to 
have a better reach. Nevertheless, a possibly stronger bound might be obtained considering 
searches in pair di-jet final states: in the RP 2 such signal can be originated by pair-produced 
(2,0) states that decay directly into SM. The analysis of this signal will be undertaken in a 
forthcoming study. 

4 Improving the LHC Sensitivity to the UED Scenario 

We display in Tab. H]the essential problem that must be confronted to improve the sensitivity 
of an LHC experiment such as CMS to the UED scenario. For all the values of ttikk studied, 
the largest cross section for producing Kaluza-Klein excitations is that for QiQi production, 
whereas the cross sections for producing a second-tier quark excitation together with a gluon 
excitation and/or another second-tier quark excitation are much smaller. Looking at the 
distributions in Figs. M and [TU1 we see that the pair production of odd states have similar 
distributions as the toy model, reflecting the similar low efficiencies for compressed spectra 
displayed in Tab. HI 

It is clear that the LHC experiments should strive to increase their efficiencies for the 
channels with the largest cross sections, particularly QiQi production. As already noted, we 
have found the sensitivities of the ATLAS m e ff search and monojet searches to be less than 
that of the CMS ot search. Therefore, the latter may provide the best starting-point for the 
optimisation of the UED search. We recall that this search was originally motivated by the 
search for supersymmetry, though the kinematical philosophy is applicable to any scenario 
with a strong hierarchy of masses and a massive, invisible dark matter particle. There may 
be scope for improving the efficiency for the UED scenario by varying the jet energy and 
rapidity analyses used in the analysis, e.g., taking into account the fact that many of the 
jets in the UED scenario are due to initial-state radiation. Improving the efficiencies for such 
channels would also benefit searches for supersymmetric models with compressed spectra. 
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H T [GeV] H T [GeV] 

Figure 11: Scatter plots for the UED model with ttikk — 700 GeV in the (H?, MHt) plane 
(left) before applying the CMS ax cut and (right) afterwards. The distribution is normalised 
to 1. 

However, a study of this possibility requires understanding the experimental capabilities in 
more detail, and lies beyond the scope of this work. 

On the other hand, we also see in Figs. |9] and [10] that the even tiers have signatures with 
very different distributions from the toy model, that are peculiar to UED models. As we 
see in Tab. H]the efficiencies for the most relevant ax and L p searches for m^K = 400,600 
and 700 GeV are relatively large for final states with (2,0) states, and are quite similar in 
all these channels for all the recurrence scales studied. These relatively high efficiencies are 
to be compared with the much lower efficiencies for the QiQi and other channels that have 
large cross sections. 

We recall that the cross sections for producing tier-2 states will be larger at 8 TeV and 
higher LHC energies than they were at 7 TeV, and that the increase will be more important 
for the even tier than for the lighter odd tier. Therefore, in designing an optimal strategy for 
the UED model it is natural to focus attention on the tier-2 states, with the aim of improving 
the reach by adapting the cuts that were optimised for searches for supersymmetric models 
that do not contain analogous states. As an aid to this effort, in Fig. [TTlwe display scatter 
plots in the (Ht, MHt) plane for the UED model with tukk = 700 GeV before applying 
the CMS ax cut (in the left panel) and after applying it (in the right panel). We see that 
the «t cut, which was largely motivated by supersymmetric models with hierarchical masses 
and corresponding MET signatures, removes large clusters of UED events with relatively 
large values of H T /MH T , notably around H T ~ 1400 GeV and MH T ~ 400 GeV and 
Ht ~ 2200 GeV and MHt ~ 100 GeV. However, a study of the cut optimisation again 
requires understanding the experimental capabilities in more detail, and lies beyond the 
scope of this work. 
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5 Summary and Conclusions 



In this paper we have explored the sensitivity to scenarios with extra dimensions of the LHC 
MET searches, which were motivated initially by searches for supersymmetry. Like the more- 
studied supersymmetric scenarios, the extra-dimensional scenarios we study have a Dark 
matter candidate, here in the form of the lightest Kaluza-Klein recurrence. However, the 
MET signal is in general reduced compared to generic supersymmetric scenarios, because the 
splitting between the masses of different states in the same tier are relatively small. We have 
used MadGraph to calculate the signals to be expected in a toy model with a single Kaluza- 
Klein quark recurrence, which we interfaced with PYTHIA and simulated for the ATLAS and 
CMS detectors using an implementation of Delphes that had been validated previously for 
supersymmetry searches. 

We then presented a more realistic model with two extra dimensions, and showed how the 
full set of first-level Kaluza-Klein recurrences would yield events with final-state leptons as 
well as MET events. We also showed that second-level Kaluza-Klein recurrences would make 
important additional contributions to the MET signals. As a result of these two additional 
signatures, the sensitivity of current LHC MET searches to the realistic model is greater 
than to the toy model, reaching into the range of masses where the lightest Kaluza-Klein 
recurrence could realistically provide all the cosmological dark matter. Specifically, in this 
first study we find that a Kaluza-Klein recurrence scale of 600 GeV is robustly excluded at the 
99.9% CL, whereas a scale of 700 GeV is disfavoured only at the 72% CL. These sensitivities 
are due principally to a combination of the CMS ax and single-lepton L p analyses, with 
opposite- and same-sign dilepton searches playing a much less important role. It will be 
interesting to extend the current study to the full LHC 8-TeV data set as soon as sufficient 
details of the MET searches are made available for calculating their sensitivity to extra- 
dimensional scenarios to be feasible. This study also showed that the reach of searches with 
MET already extends beyond the reach of dilepton resonance searches at 7 TeV, which give 
abound of 575 GeV [TBI - 

Studies of the LHC sensitivity to extra-dimensional scenarios are still much less developed 
than the corresponding supersymmetry studies, and much work remains to be done. For 
example, one potentially important class of signatures that we did not explore is pp —> X 2 
production followed by X 2 —> X±Xi, X 2 xsm decays, where X^ 2 denote generic first- and 
second-level Kaluza-Klein recurrences. The resonant production of the X 2 state takes place 
via a loop-induced vertex, and the calculation and simulation of this class of signal would 
be an interesting next step beyond the present work. We also note that there are several 
refinements of the signatures studied in this work that could be implemented, such as a 
more complete calculation of the spectra for different ratios of the sizes of the two extra 
dimensions, as well as an accurate treatment of spin effects and correlations. 

However, we think that the work presented here provides the LHC experiments with a 
benchmark that they could use in making their own interpretations of MET searches in the 
context of models with extra dimensions. In many ways, the minimal RP 2 scenario studied 
here is simpler than the CMSSM, for example, with fewer parameters. On the other hand, 
our study has demonstrated the need for a complete analysis including higher-level Kaluza- 
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Klein recurrences that have no analogue in supersymmetric models. One possible use of this 
benchmark model would be consider the design of LHC MET searches to be performed at 
higher energies after the shutdown, so that they are optimised for a broader class of models 
than popular supersymmetric scenarios. 

To paraphrase Shakespeare, "There may more things at the LHC, ATLAS and CMS, 
than are dreamt of in your current trigger, event selection and analysis philosophy" . 
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